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Abstract 

We study quasi-elastic rescattering effects in By_ (i,s DP, DP decays, wfiere P is a light 
pseudoscalar. The updated measurements of DP decays are used to extract the effective 

Wilson coefficients af^ ~ 0.90, ~ 0.23, three strong phases 6 ~ 53°, 9 ~ 18°, a ~ —88°, 
and the mixing angle r ~ 9°. This information is used to predict rates of nineteen Bg — > DP 
and Bu^d,s DP decay modes, including modes of interests in the ^/(ps program. Many decay 
rates are found to be enhanced. In particular, the Bg — > D^K^ rate is predicted to be 8 x 10~^, 
which could be measured soon. The rescattering effects on the corresponding Bu^d,s ~^ DP, DP 
amplitude ratios tb, tb^, and the relative strong phases 5b, Sb^ are studied. Although the decay 
rates are enhanced in most cases, rB,Bs values are similar to factorization expectation. 

PACS numbers: ll.30.Hv, 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



DtK- 



Color-suppressed b c decays Z}(*)°7r'' 

J3IL 

and D^K^ 5] started to emerge in 2001 (for updated measurements, see p, 7|), with 
branching ratios that are significantly larger than earlier theoretical expectations based on 
naive factorization. When combined with color-allowed B D^*\ modes in an SU(2) 
framework, the enhancement in the D^*^^tx'^ rate indicates the pres ence of non- vanishing 
strong phases, which has attracted much attention 



proposed ll| a quasi-elastic final state rescattering (FSI) picture, where the enhancement of 
color-suppressed D^h*^ modes is due to rescattering from the color-allowed D'^7c~ final state. 



This approach was also applied to study final state interaction in charmless B decays 



19|. 



The quasi-elastic approach was recently extended to i? ^ DK, DK decays [20|]. The 
color-allowed B^ D^K~ and color-suppressed B^ D^K~ decays are of interest for 
the determination of the unitary phase angle 03(or 7) = argV^*^, where V is the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix. The Gronau-London-Wyler (GLW) 21| . 



2^ and ''DK Dalitz plot' 



24| methods probe, in varying 



Atwood-Dunietz-Soni (ADS) 
ways, the interference of the two types of amplitudes in a common final state. The enhance- 
ment of color-suppressed DP modes (where P stands for a light pseudoscalar) could imply 



a larger DK rate 



15| . Since the strong interaction respects SU(3) and charge conjugation 



symmetries, the FSI in DP and DP modes should be related. It is thus of interest to study 
DP and D P modes together. 

Besides making an update with recently available data, we note that data for Bg is 
starting to emerge from the Tevetron [3| and from B factories 25|, and we anticipate more 
to come in the near future, from LHCb and other LHC experiments. Some Bg modes will 
be useful in the extraction of 7/^3 H, 28 1. It is thus timely to study Bg decays. In 
this work, we extend the scope of the quasi-elastic rescattering approach to Bg — DP, DP 



a. 



decays, as well as update our previous results using the latest Bu^ DP data 

In Sec. II we briefiy summarize and extend the quasi-elastic rescattering formula for 
Bu,d,s —>■ DP, DP decays. Numerical results are reported in Sec. III. The effective Wilson 
coefficients and rescattering parameters are obtained by using current B DP data. By 
SU(3) symmetry and charge conjugation invariance of the strong interactions, we make 
predictions on Bg DP and Bu d,s DP rates. The conclusion is then offered in Sec. IV. 
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An Appendix specifies the source amplitudes used to fit data with rescattering formahsm. 



II. FINAL STATE RESCATTERING FRAMEWORK 



We only briefiy summarize, as well as extend, the decay amplitudes obtained in the quasi- 
elastic approach for B DP, DP decays, and refer the reader to Refs. ll|] and 20| for 
more detail. 

The quasi-elastic strong rescattering amplitudes can be put in four different classes, as 
given below. For B decaying to DP with C = +1, S = 0, —1 final states, we have 
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Extending to Bg to DP decays with C = +1, S = 0, +1 final states, one has 
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For Bu,d DP decays with C — —1, S — ±1 final states, 
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And for DP decays witli C = -1, S = 0, +1 final states, 
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In these expressions, the square root of the rescattering ^'-matrix are denoted as S^^ 

(l + i7;)V2 = l + iT/, with 
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where re, jTj^ and tq are charge exchange, annihilation, singlet exchange rescattering pa- 
while fj and fj are those for DIl{8) ^ D^rji and D^i]i ^ D^rji scattering, 



rameters 



respectively [20]. SU(3) symmetry requires that %' has the same structure as %. Hence the 



T/ is basically 7^, but with rj, rj and rj replaced by r', r' anc 



r[j, respectively. We note 



that some of the above formulas were already reported in [ll|, |20|, while all formulas for 
the second and fourth cases, and some for the third case, are new. We have used charge 
conjugation invariance and SU(3) symmetry of the strong interactions, hence the r^'\ 
and •* coefficients in T^^'^ of Bu,d,s DP rescattering amplitudes are identical to those in 
of Bu,d,s DP rescattering amplitudes. 
Using SU(3) symmetry and S'^S = 1, the rescattering parameters are given by [2^ 

1 



1 



iTa = ^(3Wii-2e2^^-l), 



3 



^(ro + ^4^) = K22 - 1, (6) 



where 



COST sinr \ / e^*^ \ / cosr — sinr , 

W = W' = I , h (7) 

— sin r cos r / \ e / \ sin r cos r 

and we have set the overall phase factor (1 + iro + irg) in S to unity. This phase convention 

is equivalent to choosing the A-^o^^o^,- amplitude to be real. The r^, and in S^^'^ can 

be obtained by using the above formulas with phases (5, 9, a) reduced by half. We need 

three phases and one mixing angle to specify FSI effects in DP and DP rescattering. The 

interpretation of these phases and mixing angle in term of SU(3) decomposition can be found 



m 



20|. 



To use the FSI formulas, we need to specify A^. We use naive factorization amplitudes 
A-^ for to avoid double counting of FSI effects [ll|, [2^. And the explicit forms of A^ 
are given in Appendix A. We stress that, in our quasi-elastic rescattering approach, SU(3) 
symmetry is applied only in the DU — » DU rescattering matrix, which should hold for tub 
scale. Certain amount of SU(3) breaking effects which have to do with meson formation are 
included in the factorization amplitudes via decay constants and form factors. 



TABLE I: Branchine; ratios of various B DP and DK modes in 10 ^ units. The second column 
is ..e e.pe_.. a.a fl fl, is .a.. . .p,. The nai.e .o.e. .... 

given in third column. Fitting the experimental data shown in the first column with quasi-elastic 
FSI (fit parameters as given in Table II), we obtain the FSI fit results given in the last column. 
The factorization results are recovered by setting FSI phases in Table II to zero. 
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III. RESULTS 



In our numerical study, masses and lifetimes are taken from the Particle Data Group 
(PDG) [7|], and B to charm meson decay branching ratios are taken from We fix 

^ 0.7419. V„. . 0.2568. V. ^ 0.04166. K. ^ 0..97334. |KJ ^ 3.624 x 10- 3. and 

use the decay constants = 131 MeV, Jk = 156 MeV [7| and /dj,, = 200 (230) MeV. 

We have six parameters to describe the processes with rescattering from factorization 
amplitudes: the two effective Wilson coefficients af^ and af^ , the three rescattering phases 
S, 9 and a, and one mixing angle r in 5^/^. These parameters are fitted with rates of nine 
B decay to C = 1, 5 = 0, -1 modes, namely B D+n-, D^n- , L)°7r°, D^t], D^t]', D+R-, 
D^K' , D^K^ and D^K^ decays, given in Table I. The fitted FSI parameters are listed in 
Table II. We then use the extracted parameters to predict nineteen Bg DP (Table III) 
and Bu4,s DP (Table IV) decays. Predictions on the ratios oi B ^ D P and B —>■ DP 
amplitudes are also given. 

The errors of the fitted parameters given in Table HTl are propagated from the experimental 
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TABLE II: Fit parameters in the SU(3) FSI picture, where results are from using B D^ir , 
D+TT-, D^TT^, D^rj, D^T]', D+K-, D^R-, D+R- and D^W decay rates (Table U) as fit input. 
There is a two fold ambiguity (the overall sign of phases) in the solutions. The SU(3) phases and 
mixing angle are reexpressed in terms of the rescattering parameters r-, f-, f^. 



parameter 


result 


parameter 
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/vmin 
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K 


(0.07 ±0.01) =F (0.10±0.01)i 




i-O.lBt'oZ) ± (0.22l°:°^)z 




(0.06 ± 0.46) =F (0.97l|]:oI)i 







errors by requiring x < Ymin + '^^^ fitted values of these parameters are similar to those 
in our previous analysis ^] ^. There is a two fold ambiguity (the overall sign of the phases) 
in the solutions. We obtain x^j^/d.o.f. = 0.64 indicatin g a good fit to these modes. The 



effective Wilson coefficients af close to expectation 30, |3l|. From \r'J > \r'^\ we infer 
that exchange rescattering is dominant over annihilation rescattering. 

We show in the fourth column of Table [T] the fit output for the nine fitted B DP 
and DK modes. These fitted branching ratios (in units of 10~^) should be compared with 
data and naive factorization results given in the second and third columns. The FSI results 
reproduce the data quite well, as it should. The errors for the FSI results are from data 
only. The factorization results can be recovered by using the parameters of Table II but 
with FSI phases set to zero. Note that unitarity is implied automatically, i.e. sum of rates 
within coupled modes are unchanged by FSI. 

Our main interest here is the color-suppressed Bg decays. The predicted branching ratios 
of various Bg DP modes with C = ±1, S = 0, ±1 final states are shown in Table | 



where the second column gives naive factorization results and the third column gives the FSI 



^ We found and corrected a numerical error in our previous analysis, resulting in the value of a taking 
opposite sign. 
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TABLE III: The predictions 


on branching ratios of 


various i?s ^ modes in 


10""^ and 10"^ 


units, respectively, ine errors 


lor the 1^ bi results are 


Irom is ^ DF data only. 
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results. Again, the factorization results are recovered by using the same parameters of Table 
11 but with FSl phases set to zero, and the errors for the FSl results are from B —>■ DP data 
only. Analogous to D^n^ enhancement being fed from 5^ D^n^ rescattering, it is 

interesting to note that B^ Dfn^ with FSl rescattering to D^K'^, brings — > D^K'^ rate 
to the 10~^ level, which can be measured soon. This is helped by the absence of annihilation 
rescattering. The Bs D^r], D^r]' modes are the direct analogs of B^ i5°7r°. One can 
see that their rates are brought up to levels similar to 'B^ D^n^. Rescattering slightly 
reduces the B^ DtK~ and 5, Dt'n:^ rates. The DfK mode will be used to extract 



7/03 at LHCb |27|, while the D^r], D^rj' modes could also be useful [28 1. 

Bu,d,s DP decays are Vub suppressed, and mostly not measured yet, except Djn^. 
The quasi-elastic FSl formalism allows us to make predictions even for such modes that are 
color-suppressed. Predictions for B„,d K DP decays are shown in Table IIVI where again, 
experimental results and limits 0, [sSj are shown in the second column, and the third and 
fourth columns are naive factorization and FSl results, respectively. Same comments on 
FSl parameters apply. Agreement of the only observed B^ —>■ D^n^ mode with theoretical 
prediction is improved by including FSl effects. Bg — > DjK^ is slightly reduced, but overall, 
the redistribution of decay rates by rescattering is not very significant. 



8 



TABLE IV: Predictions for Bu^^s DP rates. Experimental results and limits j7|, |33| are shown 
in the second column, and naive factorization and ESI results are given in the third and fourth 
columns. 
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D'^K and Bs K~^, we can compare our predicted ratio 



Combining Bg 



R 



BiBs 



0.090 ±0.002, 



(8) 



with the recent experimental result of i? = 0.107 ± 0.019 ± 0.008 
agreement is reasonable. In fact, the larger rescattering of Bg -D^tt" to B 
helped enhance the ratio from the lower factorization value. 



3^ from CDF. The 
~ D^K° has 



In Table |Vl we compare our predictions for various Bg DP rates with results obtained 



in other approaches [3j, |35| that differ in the application of SU(3) symmetry. Most of our 
results agree with others. For modes with ri^'\ our results are closer to those in 3^ obtained 
using earlier data. In both approaches, U(3) symmetry is not imposed and Drji is treated as 
an independent component. Although predictions on Bg — > DP rates are similar in all three 
works, it should be note that there is a major difference between ours and the other two's 
approaches. In this work, the information obtained in Bu,d DP rescattering from data 
is used to predict not only Bg DP decays [via SU(3) symmetry], but also Bu^d,s DP 



TABLE V: 
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decays [through charge conjugation invariance of the S'-matrix]. SU(3) symmetry itself is 
not sufficient to relate -Bn,d,s DP and B^^^s DP amplitudes. Hence, in the two 
other works, which employed solely SU(3) symmetry to decay amplitudes, no prediction on 
Bu,d,s DP decays were given by analyzing Bu,d DP data. 

In Table IVll rB{Bs)^ ^b(Bs) foi^ various modes are predicted and compared with data, 
where the amplitude ratio r^^^j and the strong phase difference 5b^^-^ are defined as 
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DP) in defining Sb{b,)- Except tb^D^K'^) the 



effects from final state interaction are mild. We see that our r b{D^ K~) and 6b{D^K~) agree 
with the Dalitz analysis results of Belle and BaBar. Our r b{D^ K~) is also in agreement 
with the fit from UTjjj group obtained by using all three methods of GLW, ADS and DK 
Dalitz analysis ^3^. 



IV. CONCLUSION 

We study quasi-elastic rescattering effects in Bu^^^ g DP, DP modes. The updated 
data for nine Bu^ — > DP modes are used to extract and four rescattering parameters. 
We find the effective Wilson coefficients af^ ~ 0.90, af^ ^ 0.23, the strong phases 6 ~ 54°, 
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TABLE VI: Naive factorization and FSI results on ^b^^^^ '^Sfsi with \Vub\ = 3.67 x 10 



compared to the experimental results 
only. 



36|, 



and 



37l |. The errors for the FSI results are from DP data 
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18°, o" ~ —88° and mixing angle r ^ 9°. Since strong interaction respects SU(3) 
symmetry and charge conjugation symmetry, the formalism can be used to predict DP 
and Bu^d,s ~^ DP rates and rB{Bs)-, without refers to any experimental information on 
Bu,d,s —* DP decays, which is limited by the smallness of the corresponding decay rates. 

Our results are summarized as following: (a) The D^K^, D^rj, D^rj' rates are 

enhanced in the presence of FSI. In particular, the B^ D^K^ rate is close to 10"'^ level 
and can be measured soon, (b) The predicted Bs Dfn^ rate and the ratio of B{Bs — >■ 
DfK^)/B{Bs Z^+TT^) is in better agreement with experimental results, (c) Except the 
B^ — > D^K'^ mode, the FSI effects on rB{Bs) ^'^^ mild, (d) The predicted VBiD'^K') agree 
with data and the fit from the UTfit collaboration, while 6b{D^K~) agree with data. 
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APPENDIX A: EXPLICIT EXPRESSION OF 

As mentioned in the text, we use naive factorization amplitudes for to avoid double 
counting of FSI effects. For each final state, we have 



Ab-^D^K- = + C*/)' ^W^D+K- ^VcbV*gTf, 

4o^DOi<ro ^ ^cfeKd C"/; 4o-»D+7r- ^ ^^'^^us Ef, 

4o^DO^o = -^^^1^ Ef, 4o^D0^8 = -^^^^ (-2C/ + Ef), 

4o^z.o,, - ^ (^/ + ^/)' 42-.ntK- - (7^/ + Ef), (Al) 

where the super- and subscripts / indicate naive factorization amplitude, and 

Tf = ^af{ml-ml)fpF,^^{ml), 
Cf = ^af{ml-ml)fnFr{ml), 

Ef-^ af {ml - mDfnFr^^iml). (A2) 

pBD{BP) ^j^^ -Bu^d,s — Du,d,s{P) transition form factor and Fq~*^^ is the vacuum to DP 
(time-like) form factor. 
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TABLE VII: Form factors taken from [39|,l40|]. For 
Clebsch-Gordan coefficients are included [see Eq. (jAGj)]. 



rl form factors the mixing angle and 



Form factor 
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0.15 
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For B^, 



Du,d,sP decays, we have, 
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2c/ + e/), 



(A3) 



where, as before, the super- and subscripts / indicate naive factorization amplitude, and 



V2 



G 



72 



B) 



(A4) 



Note that we have {tf, ej) = (Tj, Ej) with D and P interchanged and Cf = Cf (without 
the interchange of D and P). 

The D^rjs and D^rji are not physical final states. The physical rj, rj' mesons are defined 
through 




(A5) 



cos d — sin "(9 
sin cos 'd 

]. Form factors are taken from [39|, |40|1, where we 
list the relevant values in Table IVIIi For —>■ r]' form factors the mixing angle and 



with the mixing angle "i? = —15.4° 
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Clebsch-Gordan coefficients are included, 



{''''D,D, 



F'^Hrnl,nJ = ("2^ + ^) Fo"^-(-?,,.J, (A6) 

where rj^ is the ss component of r] and r^' and the form factor F^'^^mjj and FQ"^"{m'^ 
are taken from [sol and [4^, respectively. 
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